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Abstract. A non-linear functional Q[u, v] is given that governs the 
loss, respectively gain, of (doubly degenerate) eigenvalues of fourth or- 
der differential operators L = d 4 + dud + v on the line. Apart from 
factorizing L as A* A + E , providing several explicit examples, and de- 
0^ ' riving various relations between u, v and eigenfunctions of L, we find u 

and v such that L is isospectral to the free operator Lq — d 4 up to one 
(multiplicity 2) eigenvalue Eq < 0. Not unexpectedly, this choice of u, 
v leads to exact solutions of the corresponding time-dependent PDE's. 
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1. Factorization of the operator L = d A + dud + v. 



^ ■ Let us assume that u and v are real- valued functions and u, v E S (R), 

^! where S (M) denotes the Schwarz class of rapidly decaying functions. Let 

^ | L be a linear fourth order selfadjoint operator in L 2 (R) 

g; (1.1) L:=d A + dud + v 

■ defined on functions from the Sobolev class i/ 4 (R). This operator is 

^ ■ bounded from below and we assume that its lowest eigenvalue E < 

is of double multiplicity and therefore there exist two orthogonal in L 2 (R) 
eigenfunctions ip + and tp^ satisfying the equation 

(1.2) Lifj = E tp. 
As shown in the appendix, the Wronskian 

(1.3) W(x) := ip + (x) ip'_(x) - ip'+(x) 

is necessarily non-vanishing, W(x) ^ 0, x E R. Let us try to factorize 

L — E as 

(1.4) A *A=(-d 2 -fd + g-f) (-d 2 + fd + g), 
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with / and g real-valued. Clearly, 

f + f 2 + 2g =-u 
9 2 -(f9 + 9j =v-E Q . 



(1.5) 



Instead of discussing these non-linear differential equations directly, let us 
express f,g,u and v in terms of the functions ip + , ip_. Straightforwardly, 
one finds that since ip + and ip- are eigenfunctions of A* A with eigenvalue 
0, we have Aip + = Atp_ = 0, which implies 



(1.6) 



fW =W 
-gW =i>'ip'L-ipiip'_=:Wi 2 , 



while [L — E Q ) ip + = (L — E ) =0 implies 
(1.7) 



uW =2W 12 -W" + e 
[v — E ) W =u W 12 + W'{ 2 - W 23 , 

where e is an integration constant and 

(1.8) ^3:=^--^- 
is expressible in terms of W and W\ 2 via 



(1.9) W W 23 = W' l2 W - W 12 W" + Wl 2 . 

Equations (11.61) say that 

W Wvi 

Since uW + W" — 2 W 12 vanishes at infinity, e has to be 0, and one finds, 
using equations dl.7t - dl.9t . that 

2 W 12 - W" 
(1.11) u=- 



W 

(1.12) v-E = T7 £+ l IJ 



W 2 V W 
Note that 

(1.13) L := AA* +E = L + 4dfd + 2fg l -ff" + f" 

will be isospectral to L, apart from E , which has been removed. To see 
why Eq is not an eigenvalue of L, let us for simplicity assume that u, v G 
Cq° (M), say that supp u, supp v C (— c, c). Then, 

ip+(x) = ax e~ KX cos (kx) + Pi e~ KX sin (kx) 
ij)-{x) = a 2 e~ KX cos (kx) + p 2 < ' ■ ; " ' ■• - 
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where E = —4 k 4 , k > 0. This implies 

W(x) = k e~ 2KX («! (3 2 ~ Pi a 2 ) 
W 12 {x) =2k 3 e~ 2KX (ai fa -ftaj) ' X > C ' 

(note that the bracket does not vanish, since ip + and ip_ are linearly inde- 
pendent.) This (and a similar investigation at the other end) implies that 

f(x) = =f2k, g(x) = —2 k 2 , for ± x > c. 

Since Lip = E ip implies A*ip = 0, we obtain 

ip" -2kV' + 2k 2 ^ = 0, x > c. 

It clearly follows that ^ cannot be in L 2 (R) unless it vanishes identically. 

Before giving some explicit examples, let us make some comments con- 
cerning the problem of actually finding / and g, or ip + and ip-, when u and 
v are given. Instead of solving the non-linear system (11.51) . or the spectral 
problem (11.21) . one may also try to solve the Hirota-type equation which 
follows from (TTTTTT) . dTTT2b 

(1.14) i{v -E )=( w + uj + 2 (— + + , 

and which for u = reads 

4 (t> - £ ) W^ 2 = 2 (W" W - W" W") + W" 2 . 

Once 0) is obtained, / and g can be given by the equations (11.101) . 

With / and g defined in this way O, equation (11.51) is satisfied and the 
factorization (11.41 ) is valid. 

Note also the following: the functions ip + and ip- are solutions of A ip = 

0, i.e. 

-ip" + fiP' + gip = o. 

By writing 

il>± = Vw</> ± , 

one finds that <fr + <f>'_ — (p' + <p- = 1 and that <p± are (oscillating) solutions of 
the equation in Liouville form 

,.. / 3 fW'V 1 W"\ , 

1. e. associated to a second order self-adjoint diffential operator. 
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2. Addition and removal of eigenvalues. 

Although adding and removing eigenvalues may be thought to be a pro- 
cedure that can be read both ways (symmetrically), the steps involved are 
actually quite different in both cases (in particular, it is not yet clear, which 
conditions on u and v allow for the addition of a doubly degenerate eigen- 
value below the spectrum of d 4 + dud + v). Let us therefore 'summarize' 
them separately, in both cases starting from a given operator 

L n := <9* + d u n d + v n , n E N, 

and the equation (11.141) with u, v replaced by u n , v n . This equation shall be 
referred to as (11.141) ,,. 
Removal of eigenvalues: 

1. Solve (TLUL (with E -> E& n) = -4k*) for W n := W 

(-> at infinity) and define W$ as \ (W n u n + W"), as is natural in ac- 
cordance with equation (11.1 1L . Alternatively, if ip^ are known, calculate 
W n and W[f via their definitions, i.e. as 

w§> =^ ) > w,/ -^ )/ V w '. 

2. Define /„ and g n according to dl.lOL . thus solving the system (I1.5I) . and 
obtaining the factorization 

L n = A n A n — 4k*. 

3. The operator 

L n — A n A n — 4K n =: L n -i 

will then be isospectral to L n apart form the lowest eigenvalue = —4k* 
(of multiplicity 2), which has been removed. 
Addition of eigenvalues: 

1. Solve dlT4L (with£ -> E^ n+1) = -4k* +1 ) for W n+1 := W ~ 

e ±2Kn+lx , as x — > ±oo, i.e. W n +i diverging at infinity and non-vanishing for 
finite x. (As mentioned above, conditions on u n , v n ensuring the existence 
of W n+ i are still unclear.) 

2. Define W n +i ■= p^r - > which will then solve the (more complicated 
looking) equation 

W' A W"" W"'W' W" 2 W"W' 2 

(2A) 40 T^- 2 ir + 14 — + 13 t^- 64 — + 

W' fw 12 (W'\ » 
2h" + ii 2 -2m'— + 2u( — -2( — | ) = Hi,. 1 ., 
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(with u,v^ u n , v n and k — > n n +i)- hi fact, (12.11) is equivalent to 

- 2f + Qff" + 7f 2 - 8ff + f + 2u (f - 2f>) - 2u'f + u 2 + 2u" 
= 4v + 16k 4 

(via / = =: /„ + i, m, v — > u n , v n and k —>■ K n+ i j that arises in the 

factorization of L n+ i. 

3. Write 

L n = A n+ iA* n+1 — 4k 4 1+1 

(implying 2 g n+1 := 3 - f 2 +1 - u n ). 

4. Then, 

L n +l '■= ^-n+l^-n+l — 

will be isospectral to L n apart from having one additional (doubly degener- 
ate) eigenvalue E^ 1 ^ below the spectrum of L n . 

3. A NON-LINEAR FUNCTIONAL Q AND A SYSTEM OF PDE'S 
ASSOCIATED WITH THE OPERATOR L. 

As observed 100 years ago [4|, the operator L = <9 4 + dud + v has a 
unique 4'th root in the form L 1//4 := d + Y^T=i h( x )d~ k . Define M to be 
the positive (differential operator) part of any integer power of L 1//4 . Then 
it is well known, that 

L t = [L,M], 

where L t is the operator defined by L t ip = du t dtp + v t (f, consistently 
defines evolution equations (for u = u(x, t), v = v(x, t)) that have infinitely 
many conserved quantities (i.e. functionals of u and v, and their spatial 
derivatives, that do not depend on t). We shall make use of this by letting 

M :=8(L 3/4 ) + = 8d 3 + Qud + 3u, 

and focusing on the quantity 
(3.1) 

Q[u,v) := j ^48v 2 + ^u 4 -12n 2 w-40Mv"-13u'u ,2 + 9w" 2 ^ dx. 

This quantity does not change when u and v evolve according to 

u t = 10 u"' + Quu' -24V 
v t =3 (u'"" + u v!" + v! u") - 8 v'" -6uv'. 

Formula (11.131) for L = d 4 + d u d + v implies that 

u — It = 4 /' 

v-v =2fg'-ff" + f". 



(3.2) 



(3.3) 
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By using the asymptotic properties of / and g (f — > q=2 k, g — > —2k 2 , as 
x — > ±00), one can show that 

2 9 

(3.4) 5Q := Q[u,v] - Q[u,v] = -32tc 7 y. 

^making SQ = —2 (4K 4 ) 7 ^ 4 j. This result is similar to that for 

Schrodinger operators [ 1 1 and reflects the loss of a doubly degenerate eigen- 
value E = —4 k 4 , when going from L to L. The constant in the right hand 
side of (13.41) is related to the semiclassical constant appearing in the trace 
formula for a fourth order differential operator considered in Q. 
The proof of (13.41) . just as the derivation of d3.ll) . involves very lengthy 
calculations. When deriving (13.41) one uses (11.51) and (13.31) to write the ex- 
pression for SQ as an integral of terms involving only the functions / and 
g, and their spatial derivatives. The crucial step is to note that the integrand 
is a pure derivative of x, i.e. SQ = J Q ' dx for some function Q, which 
makes it possible to evaluate the integral solely from the limits of / and g 
at infinity. Thus, to compute SQ, one selects the terms in Q which are free 
of derivatives, as those are the only ones that contibute. The terms in Q still 
containing derivatives, for instance the ones quadratic in g and linear in /, 

y T(96 - 48) g 2 f" - 2 • 96 / g' g" - 8 • 12 g" f ■ 2g - 4 • 40 /'" g 2 

+ 160 / g'f - 16 • 26 /" gf g - 16 • 13 g' 2 f'^j dx, 
give zero. 

4. Some examples. 
Example 1. The operator 

L = d A - 5 d 2 + d =- d V- = A * A - 4 

cosh x cosh x 

with 

A = -d 2 - 3 tanhxd- 2 
has 2 linearly independent eigenfunctions with eigenvalue E = —4, 

sinh x 



cosh 2 x ' cosh 2 x 

One can easily check that A ip± = and that it, v are reflectionless, as 

Z = aa* - 4 = a 4 - 5 a 2 
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(note that ip + and have different fall-off behaviour at oo and that 

W(x) = cosh~ 3 x). 
Example 2. The operator 

1 40 88 

L = a 4 + 16 d 5- d + ^- = A* A - 64 

cosh x cosh x cosh x 

with 

A = -5 2 - 4 tanh x 9 - ? 



cosh 2 x 

has 2 linearly independent eigenf unctions with eigenvalue i? = —64, 

. , x cos 2x , , . sin2x 

M x ) = — T2-> V-w = — p-- 

cosh x cosh x 

One easily verifies that A = 0, and that 

40 



L = AA* - 4 = d 



i 



cosh 2 x 

A computation gives that 

oio _ 2 21 / 9 9 

Q = _.4-687, Q = 2 10 -100, = — — ( = -32 (k = 2) 7 y 

Example 3. The operator 

L = a 4 + (45 ^ 4 - 40 ^ 2 ) = A*A - 4 

with 

W = ^ 2 := W 12 = 2^ 2 -3^ 4 

cosh x 

and 

A = -d 2 -2 tanhx<9-2 + 3^ 2 
has a doubly degenerate eigenvalue i?o = —4. One easily verifies, that 

L = <9 4 - 8 d ^ 2 a + 25 ^ 4 - 16 ^ 2 . 

Example 4. The operator 

L = d 4 -d 2 + 4c>— J^d+^ \^- = A*A 

cosh x cosh x cosh x 

with 

A = — d 2 — tanh x9 ~— = d (—d — tanh x) 

cosh x 

has a unique ground-state £o = with eigenfunction 

if)(x) = — l — . 
coshx 

The second solution of Aip = is ip = tanhx ^ L 2 (R). One easily 
verifies, that 

Z = d 4 - d 2 . 
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Example 5. For any k > 0, the operator 

L = d A + dud + v 

with 



u{x) =2(l + f)^(|) 



where 



cosh x ' 

has a doubly degenerate ground-state, E = —4, with eigenfunctions 

/ i \ k 

iP ( ±\x) = e ±lx 



1 



cosh | 



5. FOLLYTONS. 



In order to find u and t> such that L = A* A + E is 'conjugate' to the free 
operator L = d 4 =: L one has to solve (11.51) with w = v — 0. Eliminating 
g and writing £? = — 4k 4 one obtains the ODE 

2f"' + 6ff" + 7f' 2 + 8f'f 2 + f 4 = 16^. 

One may reduce the order by taking / as the independent variable, and 
F{f) '■= f as me dependent one, yielding 

2 [F" F 2 + F' 2 F) + QFF' f + 7 F 2 + 8F f + f = 16 k 4 , 

but both forms seem(ed) to be too difficult to solve. By using (11.141) . how- 
ever, it takes the form 

16k 4 W 2 = 2 {W"" W - W" W) + W" 2 , 

in which it is easier to see the solution Q 

W = const ■ (V2 + cosh (2 « x) \ . 

Correspondingly, 

p W' sinh (2 k x) 



W y/2 + cosh (2 « 



x 



As interchanging A* and A (as far as / is concerned) only changes the sign 
of/, 

_ . . sinh (2 k x) 



(5.1) 



v K 



(5.2) 
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The Wronskian of the two ground-states ip± (of L = d i + dud + v, conju- 
gate to L = <9 4 ) is simply the inverse of W, i.e. (choosing the constant in 
W to be 1), 

W(x) = —= =: x (2 k x) . 

V2 + cosh (2 k x) 

The function g is given by 

g=~(3f'-f 2 ) = -2K 2 (1 + V2W-2W 2 ). 
Insertion into equation (11.51) yields the reflectionless 'potentials' 

= 16k 2 (V2W- W 2 ) 

-- 16 k 4 (V2W -12W 2 + 16V2W 3 -8W 4 ) 

with L = d A + du K d + v K having exactly one doubly degenerate negative 
eigenvalue —4 k 4 . While in most other examples we scaled k to be equal 
to 1 it is, in this case (due to the appearence of 2k in W) easiest to choose 
k = i, i.e. to take 

u =4( v / 2x-X 2 ) 
v = (V2 X -12 X 2 + 16V2X 3 -8X*) 
and, when needed, use formulas like 

X " =X (l-3V2 X + 2 X 2 ) 

X 12 =X 2 (l-2V2 X + X 2 ) 

X'" =x' (1-6V2X + 6X 2 ) 

X"" =X (l-15v / 2x + 80 X 2 - 60v / 2x 3 + 24 X 4 ). 

(Note that redefining x by a factor of —\pl would make all the coefficients 
positive (integers)). These formulas are useful when checking that n(x+4 1) 
and v(x + 4 1), with u and v given by (I5.2I) . are exact solutions of the non- 
linear system of PDE's (13.21) (just as u K (x + 16k 2 t), v K (x + 16k 2 t)). 

Appendix A. W ^ 0. 

We shall prove here that the Wronskian type function defined in (11.31) 
never equals zero. 

Theorem A.l. Let if)± be two orthonormal eigenfunctions of the operator 
(11.11) corresponding to the lowest eigenvalue E of double multiplicity. Then 

(A.l) W[i(j + ,ip-](x) := ip+(x) ip'_(x) - ip-(x) ip' + (x) ^ 0, xeR. 
In order to prove this result we need a simple auxiliary statement. 
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Lemma A.l. Let E be the lowest eigenvalue of the operator L and let tp G 
L 2 {R) be a solution of the equation (11.21 ) satisfying ip(x ) = ip (x ) = 
for some Xq G R. Then ip(x) = 0. 

Proof. Indeed, the function 

~ \i/>(x), ifx< x , 

ip(x) — < 

I —ip(x), if x > x , 

is linear independent with tp. Since ip(x ) = ip (x Q ) = we obtain 
(Lip,ip) = E \\ip\\ 2 . Then ip is also an eigenf unction of the operator L 
with the eigenvalue E . Consider now the linear combination 

ipx(x) = -tp" (x )ip(x) - xp" (xq)tP(x). 

Obviously ipi( x o) — ipi{%o) — ^\"{ x o) = 0, ip\ G L 2 (R) and tpi satisfies 
the fourth order differential equation Lip± = E ipi. Being overdetermined, 
■01 = which also implies ip = 0. □ 

Remark. In Lemma lATTl the conditions ip(x Q ) = ip' (x ) = split the prob- 
lem for the operator L in L 2 (R) into two Dirichlet boundary value problems 
on semiaxes L 2 ((x , oo)) and L 2 ((— oo, x )). Therefore, the lowest eigen- 
value moves up. 

Proof of Theorem \A.1\ 

a. Let ip±, be two orthonormal eigenfunctions corresponding to the lowest 
eigenvalue E of the operator L. The functions ip + and ip- cannot vanish at 
the same point. Indeed, assume that they do. Then there is a point x such 
that ip + (x ) = ip-(x ) = 0. If in addition we assume that say ip + (x ) = 0, 
then by Lemma |A~T1 ^ + = and we obtain a contradiction. Therefore we 
can assume that ip ± (xo) ^ 0. Introduce a new function 

2 (x) = ip'_(x )ip + (x) - ip + (x )ip-(x). 

It is a non-trivial eigenfunction of the equation (11.21) satisfying ^2(^0) = 
ip 2 (x ) = 0. By using Lemma lATTl again we find that tp 2 = which cannot 
be true because ip + and ip- are linear independent. 

b. Consider now the following pair of complex functions 

m ± ( x ) = -,p + {x) ± iip„{x) =: ip(x)e ±i ' t ' {x) . 
By using a. we observe that ip never vanishes, ip(x) 7^ 0, x G R. Besides 

= (ip + + iip-){ip+ - - (ip+ + iip-)'(^+ ~ iip-) 
= -2iW[ip+,ip_] = -2i(P'ip 2 . 
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Thus, in order to prove Theorem IA. II it remains to prove that <f> ^ 0. As- 
sume that there is x such that (x ) = and consider 

Clearly $(a?o) = $ (xo) = and by using Lemma |A~T1 we obtain $ = 
which contradicts the linear independency of the functions 
The proof is complete. □ 
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